Directional emission of stadium-shaped micro-lasers 
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The far-field emission of two dimensional (2D) stadium-shaped dielectric cavities is investigated. 
Micro-lasers with such shape present a highly directional emission. We provide experimental evi- 
dence of the dependance of the emission directionality on the shape of the stadium, in good agree- 
ment with ray numerical simulations. We develop a simple geometrical optics model which permits 
to explain analytically main observed features. Wave numerical calculations confirm the results. 
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The field of quantum chaos has widely broadened over 
the last two decades 1]. Generally, it relates the quantum 
behavior of a broad diversity of systems with their classi- 
cal features. In this context, lossless billiards are models 
of great interest, providing a broad variety of dynamical 
systems by changing the shape of the boundary. More- 
over they are accessible to experimental studies 0, 1^- 
In the same time, the influence of loss or noise on open 
quantum systems is also being investigated 0, with em- 
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phasis on the behavior of quantum systems with 

chaotic classical dynamics 5,"^, V]. Flat micro-lasers ex- 
hibiting different boundary shapes are relevant examples 
of open billiards with a coherent output coupling 

Here we focus on polymer micro-lasers of the Buni- 
movich stadium shape. Such billiard is the archetype of 
chaotic systems ,2j . Stadium billiards look like a rectan- 
gle of length 21 between two half-circles of diameter 2r, 
accounted the form ratio l/r (see Fig. 0. Though all pe- 
riodic orbits are unstable, dielectric micro-cavities with 
this shape are well-behaved lasers with highly directional 
emission in the far-field pattern Aside from quan- 

tum chaos studies, such micro-resonators are also being 
considered for applications in integrated optics and bio- 
logical sensors pTj . 

This letter is devoted to the investigation of the direc- 
tional emission versus the form ratio as a way to infer 
general behavior for such dynamical systems. First we 
present experimental results and show a good agreement 
with ray numerical simulations. Then we develop a sim- 
ple ray model which permits to describe analytically main 
observed features. All these results are finally confirmed 
by electromagnetic numerical calculations. 

The micro-laser cavities are etched in a thin layer 
of a passive polymer matrix (PMMA: polymethyl- 
methacrylate), doped with a guest laser dye (DCM: 
4 — dicyanomethylene — 2 — methyl — 6 — (4 — dimethyl- 
aminostyryl) — AH — pyran). Our versatile fabrication 



FIG. 1: Left: Optical microscope image of a stadium-shaped 
micro-laser with l/r = 1 and r — 30 jj,m. Right: Nota- 
tions. I is the half-length of the rectangle and r is the radius 
of the half-circles, s is the curvilinear coordinate along the 
boundary, x is the incident angle and 6 is the outgoing angle 
measured from the main axis of the stadium. 



process ensures for broad variations in size {r — 10 /im 
to 50 /im), shape {l/r = to l/r — 5) and thickness (be- 
tween 0.4 /im and 0.7 /im to avoid vertical multimode be- 
havior) with quality factor Q greater than 6000 10]. The 
microlasers are pumped uniformly from above at 532 nm 
with a frequency doubled picosecond Nd:YAG laser and 
emit from their sides with A ~ 600 nm. Such 2D dielec- 
tric billiards present an effective refractive index n of 1.5 
and are operating in the semi-classical regime (r/X vary- 
ing from typically 20 up to 100). The emitted light is 
collected in the far-field and coupled to a multi-channel 
detector via a spectrometer. A typical experimental spec- 
trum is shown in Fig. |21 (bottom). The sample is fixed 
on a rotating mount and the intensity emitted in a given 
direction is obtained by summing over all pixels of the 
spectrum. Far-field spectra are performed over 360° by 
5° steps with the aperture angle of about 15°. 

The intensity of the laser emission versus the angle 
(see Fig. ^ (right) for notations) is displayed in Fig. |3 
(top). The emission exhibits a high directionality along 
four directions, symmetrical with respect to the 0° and 
90° axis, according to the obvious symmetries of the sta- 
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FIG. 2: Top: Experimental far-field emission in the plane of 
a 2D stadium-shaped micro-laser as a function of the polar 
angle for l/r = 0.5 (solid line) and l/r — 1 (dashed line). 
Bottom: Experimental far-field spectrum for l/r = 2 and 
r — 17.5 nm. 

dium shape. The direction of maximal emission depends 
on the l/r form ratio. Experimental results are summa- 
rized in Fig. Olfor form ratios ranging from l/r = 0.5 to 
l/r — 3. They represent an optimum for reproducibility 
according to our 5° precision interval. 

Even if stadium-shaped billiards are fully chaotic sys- 
tems, these micro-lasers present experimentally a highly 
directional emission in the far-field pattern. The direc- 
tion of maximal emission depends on the l/r form ratio 
and a variation of more than 30° has been measured be- 
tween l/r — 0.5 and l/r = 3 (Fig. I^Jwhat can be of great 
interest for applications in integrated optics. 

To account for these experimental results, we first per- 
form usual ray numerical simulations consisting in prob- 
ing a large number of randomly defined rays (^ 10^) 
with starting points and initial directions uniformly dis- 
tributed over the whole phase space. Each ray propagates 
along a straight line and is totally reflected at the bound- 
ary until the incident angle x becomes smaller than the 
critical angle, Xc = arcsin(l/n), thus allowing the ray 
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FIG. 3: Direction of maximal emission Q versus form ratio 
l/r. Squares represent experimental results. Solid line in- 
dicates dominant contribution of ray numerical simulations 
and dotted line is the position of the largest satellite peak. 
Dashed lines are lens model predictions with m = (lower 
line) and m = — 1 (upper line). Inset: Comparison of ray 
(solid line) and wave numerical simulations (circles). Dashed 
lines are lens model predictions with (from bottom to top) 
m = 0,-1,1,-2. 

to escape by refraction 15] . The stadium-shaped bil- 
liard is fully chaotic, so almost each ray escapes after 
propagating over a finite distance. Typical histograms 
of the outgoing angles are plotted in Fig. ^ To remove 
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FIG. 4: Histograms of outgoing angles obtained by ray nu- 
merical simulations with threshold equal twice the billiard 
perimeter for l/r = 0.5 (solid line), l/r = 1 (dotted line), and 
l/r = 2.5 (dashed line). Pictures normalized to unit maxi- 
mum. Inset: The same for l/r — 0.5 but without cut off. 

the background, only those rays which propagate over 
more than a finite threshold distance (twice the billiard 
perimeter in Fig. 2)) are taken into account. Histograms 
are indented by peaks which reveal the structure of un- 
derlying unstable manifolds 0, An histogram is 
calculated for each form ratio value and the position of 
its maximum is plotted versus l/r in Fig. 13 (solid lines). 
The plot is in excellent agreement with real experimental 
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results as seen from Fig. |3| 

For some l/r values ray simulations reveal the exis- 
tence of a few peaks with comparable amplitudes (cf. 
Fig. 0)). Dotted line branch in Fig. 13 indicates the 
most prominent satellite peak. Nevertheless dominant 
experiment contribution is assumed to correspond to the 
upper branch because its associated intensity (integral 
below the peak) is larger than the other one (cf. solid 
line curve in Fig.QJ. 

To explain main features of experiments and ray sim- 
ulations, we develop a simple model based on geometri- 
cal optics. To minimize refraction loss, rays propagating 
near the boundary are privileged. So we focus on rays 
which, after being totally reflected on the left half-circle, 
come to bounce on the right one (see Fig. 0). Long-lived 
trajectories can be characterized by their density on the 
Poincare surface of section with coordinates s and sinx. 
Our main assumption is that this density is uniform in 
the allowed part of phase-space 



P(s,sinx) 



|.s|)e(|sinx|-l/7i) (1) 



where Q{x) is the Heaviside step function. 

Long-lived trajectories in open chaotic systems are 
concentrated near unstable manifolds of confined peri- 
odic orbits ^ J2i . So the true density is a fractal whose 
correct determination requires the knowledge of large 
number of periodic orbits. However, finite experimen- 
tal resolution and quantum effects will smooth out small 
details and as the first approximation the assumption 
looks quite natural. Below we show that it leads to sim- 
ple analytical formulas in a good agreement with direct 
ray simulations. 

Let us consider point sources on the left half-circle 
emitting light according to Among these, some rays 
escape by refraction on the right half-circle and one can 
deduce the histogram of outgoing angles. The depen- 
dence of the histogram maximum on I jr reproduces well 
experiments and ray simulations, validating the assump- 
tions. As the resulting curves are closed to lens model 
predictions (see below) we did not present them. 

To further simplify this approach, we consider the right 
half-circle as an ideal lens. In this approximation all rays 
emitted from a point source are focused exactly as rays 
close to the line connecting the source and the center 
of the right half-circle (i.e. spherical aberrations are ig- 
nored) . So the direction of emission d corresponds mainly 
to this line and can be calculated analytically from purely 
geometrical considerations. It depends on the emission 
angle x ^"^^ on the number |m| of rebounds on straight 
boundaries (m > - resp. m < - means first bounce 
on the high - resp. low - straight boundary): 



{llrf 



arcsm 



smx 



This line exists provided x\ < sinx < X2 where xi = 1/n 
and X2 = 21/ y^AP + (2m + l)^r^. Formally the channel 
with given m is open when xi < X2 or l/r > |2m -I- 
l\/2\/n'^ — 1 but it may exist for smaller l/r as well. For 
a given l/r, the mean emission angle, 9, is deduced by 
averaging over the above interval 
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X2 ~ Xi 



6l(x)dsinx 



(2) 



Within the lens model the emission intensity is propor- 
tional to the angle with which the right half-circle is 
visible from a source point. Direct geometrical calcula- 
tions show that the dominant contribution for l/r < 3.5 
corresponds to m = channel. Channel with m — —I 
has smaller but comparable intensity. Considering 
experimental precision, other channels can be ignored 
but they are visible in ray simulations (cf. a small peak 
around 30° in the histogram for l/r — 2.5 in Fig. 0J. 
The directional emission predictions for channels with 
m — and m = —1 are indicated by dashed lines in 
Fig. The excellent agreement between these curves, 
experiments and direct ray simulations confirms that the 
assumption is relevant to interpret analytically the 
emission of the stadium-shaped cavity in the geometrical 
optics domain. 

For completeness we perform also numerical calcula- 
tions in the wave domain. In the simplest approxima- 
tion, the electromagnetic field is evolving in a passive 
cavity with a defined polarization. In this case, the elec- 
tric field (TM polarization) or its magnetic counterpart 
(TE polarization) is represented by a scalar wave func- 
tion -0 obeying the two dimensional Helmholtz equation 
[(A + n^gfc^) ip — where refractive index rii = n in- 
side and rie = 1 outside the cavity] with specific bound- 
ary conditions (see e.g. ^11,]). To find eigenvalues {km} 
and eigenfunctions {ipm} of quasi-stationar y s tates we 
use the standard boundary element method 14^. In this 
approach the wave functions inside and outside the cav- 
ity il'i,e are expressed in terms of single layer potentials 



<j> ds /ii,e(s) H^^\kr 



.\r-rs\) (3) 



s is the curvilinear coordinate along the boundary. Here 
Hq^^ is the Hankel function of the first kind. The choice 
of this Hankel function outside the cavity is dictated by 
outgoing conditions to infinity 



ikr 



m 



(4) 



2Vm2 + {l/rf 



and inside the cavity by stability of numerical algorithm. 

Following natural filtering by the lasing effect, only 
modes with small losses are considered, i.e. modes with 
an imaginary part of the wave-number closed to zero. 
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Most of the well-confined wave-functions present a clear 
directional emission in the far-field pattern, even for low 
\kr\ (see Fig. |31 (right)). Then it is possible to find po- 
sitions of directional dominant peaks for each l/r and 
different eigenvalues which are plotted in Fig. |3| (inset). 
Each point at fixed l/r corresponds to the averaging of 
maximal emission angles over a representative set of well- 
confined eigenfunctions. The error bars corresponds to 
the mean standard deviation. As was predicted, dom- 
inant emission direction corresponds to m = and 
m = —I channels of the lens model. It is interesting 
to notice that channels with m — I and m — —2 are also 
well visible. 

So good agreement between geometrical and wave op- 
tics does not usually appear for closed chaotic systems. 
The emission directions of these quantum dielectric bil- 
liards behave largely as if they were classical ones, gov- 
erned by geometrical optics rules. However there is no 
loss of coherence as in open systems generally studied in 
the field of decoherence |j| : experimental spectra present 
clear peaks (see fig. 12 (bottom)) and no noise is intro- 
duced to couple light with the outside of the cavity. Ac- 
tually the typical distance L/r covered by a photon - 
of the order of 2Q/n\k\r - ranges between 16 and 80 in 
our experiments. This magnitude suggests that inter- 
ference effects should be taken into account whereas we 
have demonstrated that a geometrical optics approach is 
sufiicient to predict the overall directional emission. The 
lens model shows that it is mostly connected with the 
existence of only a small number of well defined open 
emission channels. Under higher precision calculations 
and experiments, each channel should split into many 
narrow peaks and quantum phenomena should be more 
important. 
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FIG. 5: Left: Imaginary part of quasi-stationary eigenval- 
ues versus their real part. These eigenvalues correspond 
to solutions of the Helmholtz equation with TM boundary 
conditions for dielectric stadium with l/r — 2 and refrac- 
tive index n=1.5 antisymmetric with respect to both axes. 
Right: Far-field emission pattern of eigenfunctions: |/(^)P 
in @ normalized to unit maximum. From top to bottom 
{l/r^O.3, i?e(fer)=39.31) ; {l/r^l, i?e(fcr) = 13.88) ; {l/r=2, 
i?e(fcr)=23.29). 

In summary, we present experimental results about 
directional emission of stadium-shaped micro-lasers and 
demonstrate that they are well described by simple ge- 



ometrical optics approaches. We develop a lens model 
which predicts analytically main features of directional 
emission. Moreover we point out that well-confined wave- 
functions show a good agreement with these classical pre- 
dictions while coherent properties are not destroyed. 

These kinds of predictions seem to be applicable to a 
more general class of chaotic dielectric billiards. Experi- 
ments and calculations are thus under progress to explore 
other micro-cavity shapes like cut disk and cardioid. 
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S. Brasselet for experimental and technological support 
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